CONTINUOUS VARIABLE TOMOGRAPHIC MEASUREMENTS 
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Abstract. Using a recent result of Albini et al. [T to represent quantum homodyne tomog- 
raphy in terms of a single observable (as a normalized positive operator measure) we construct 
a generalized Markov kernel which transforms (the measurement outcome statistics of) this 
observable into (the measurement outcome statistics of) a covariant phase space observable. 
We also consider the inverse question. Finally, we add some remarks on the quantum theoreti- 
cal justification of the experimental implementations of these observables in terms of balanced 
\ homodyne and 8-port detection techniques, respectively. 
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There are two main approaches to continuous variable quantum tomography. The first one 
c3 ■ apphes measurements of rotated quadratures Qg, 9 G [0, 27r), the second one uses a measurement 
D^. of the phase space observable G|o) generated by the vacuum state |0). In both cases the state 
p of the measured system can actually be reconstructed from the corresponding measurement 
statistics, either from the quadrature distributions p^'',6' G [0,27r), or from the phase space 
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1. Introduction 
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A direct comparison of the two tomographic methods has suffered from the fact that in the 
first case one needs to measure, in principle, infinitely many different, mutually noncommuting 
observables Q^, 9 G [0, 27r), whereas in the latter approach a measurement of a single observable 
is sufficient. For each state p the random sampling of the distributions p^", 9 G [0,27r), 
can, however, easily be amalgamated into a probability bimeasure, which, when taken all 
together, define a unique informationally complete positive operator measure Eht, the POM of 
the quantum homodyne tomography [I]. 

In this letter we construct a generalized Markov kernel which turns (the statistics of) the ho- 
modyne tomography observable Eht into (the statistics of) the phase space (Husimi) observable 
G|o). We also demonstrate that there is no simple inverse construction. Finally, we briefly com- 
ment on the theoretical background of the physical implementation of the observables Eht and 

G|o) via balanced homodyne detection and 8-port homodyne detection tehniques, respectively. 
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2. Notations and the two POMs 

Let hni n G N, denote the (normalized) Hermite functions in 7i = L^(]R), the corresponding 
number operator, and let Q and P stand for the usual position and momentum operators, 
respectively. Let Q denote the spectral measure of Q so that the quadrature observable Qe is the 
spectral measure of the rotated position operator e^^^ Qe~^^^ . For any state p, positive trace one 
operator on 7Y, we let p^" denote the density of the probability distribution X i-^ tr [pQe^X]]^ 
that is, tr[pQ0(X)] = /^pQ»(x) dx. 

For each state p, the random sampling of the distributions p^" can, indeed, be amalgamated 
into a probability bimeasure i3([0,27r)) x BiW) 3 (e,X) ^ /^tr [pQe(X)]£ G [0, 1], which can 
uniquely be extended to a probability measure on i3([0, 2ti) x M) (we use the notation B{- ■ ■ ) 
for the Borel sets). As shown in [IJ, these probability measures define a unique normalized 
positive operator measure, the POM Eht, with the propertyl3 

(1) Eht(e X X) := ^ / Qe{X)de = ^ [ e^^^Q(X)e-^^^rf^, 6 G S([0,27r)), X G B{R). 

For any state p, we let p^^^ denote the density of the probability measure Z \—>- ti [pEht(^)] 
with respect to ^d9dx, so that p^^^{6,x) = p^*(x). Since the Radon transform of the Wigner 
function Wp of p gives the quadrature distributions, that is (7lWp){9,x) = p^''(x), for almost 
all 6 G [0, 27r), X G M, one thus haC 
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(2) tr[pEht(exX)] =— / / p^^^^i9,x)dedx = — / 
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Let Wqp = e^e~*'^^e*^'5 denote the unitary Weyl operators. As well known, see e.g. [71 [211 111 
[TT] . any normalized POM G which satisfies the covariance condition with respect to the Weyl 
representation of the phase space translations, has an operator density {q,p) ^ -^WqpKW*^ 
defined by a unique positive trace one operator K, that is, for any state p, 

(3) tr [pG(Z)] = ^ ^ tr [pWqpKW*J dqdq. 

^ Obviously, one could equally well define the POM Eht(0 x X) := i Qg{X)d6 for phase variable in [0, tt). 
Then Eht(exX) = iEin (On [0, tt) x X) + iEht([(en [tt, 27r)) - tt] x {~X)) for all 6 G 6([0,27r)) and X e B{R). 
To have a direct use of some of the results of [I] we choose to work with Eht, which, however, contains some 
redundancy in giving a double covering of . 

^The third equality is obtained whenever p is such that Wp G L^(]R^), see, e.g. [T]. 



In particular, if K = \ho){ho\, then the phase space observable Gk generated by K is denoted 
by G|o), and the density p^i"> of the probability measure Z t— > tr [pG|o)(2')] is just the Husimi 
Q-function of the state p, that is, p^^°^{q,p) = {Wgpho \pWqpho), where Wgpho = \z) is the 
coherent state, with z = -^{q + ip)- 

3. Construction of the kernel 

We shall construct a generalized Markov kernel M|o) (defined on M? x ([0, 27r) x M)) such that 
for any state p, the density p^i"> is obtained from the density p^i^* as follows: 

(4) p^i°Kg,p) = Y r I ^ror(^.^)p'"(^>^)^^^^' 

or, equivalently, for the ROMs, 

G|o)(Z)= r II M^^{d,x)dqdj\d£^,{d,x)= r [ M\o){Z-e,x)dEu{9,x). 
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where Z C is compact (so that we can change the order of integration). 

In a forthcoming paper [TH] a general method is developed to construct such kernels for 
arbitrary phase space observables Gk- This general method shows that the kernel can be 
obtained as a derivative of the Hilbert transform of the density of the generating operator K. 
In the special case of K = \ho){ho\ one may confirm this fact by a direct computation. We shall 
do so next. 

Before doing that we note, however, that essentially the same computations have been carried 
out in a different context of quantum-state sampling where the matrix elements of the state, the 
density matrix, are constructed form the quadrature distributions in terms of a kernel function. 
In that connection these functions were called the sampling functions or the pattern functions, 
see, e.g. [151 Sect. 5] and the references therein. 

The Hilbert transform of the function /iq, /lo(a^) = can easily be computed and one 

observes that it is the Dawson integral 

(5) daw(x) = e-"' /" e''dt. 
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Following the method of pS] we now define the function 

d 

(6) M-'^P{9,x) = 2— daw(x - gcos^ -psin^). 



which is known [10] to be a bounded analytic function vanishing at infinity and having the 
values 
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M^fie, ^) = Z2 2fc(4)! ^'' {x-qcose-p sin ( 



It can easily be seen that 2 ddaw{x)/dx is not a positive function (see for instance the picture 
in page 114 of [13]) so that J^\o){Z;9,x) is not a true conditional (or transition) probability; 
therefore we call A4\o){Z; 6,x) a generalized Markov kernel. 
We go on to show that (for compact sets Z) 



(8) 



G|o>(^) 



2lT 



— / M, 

2n 



x)dqdp 



dEu{9,x) 



It suffices to show that the coherent state expectations of both sides of this equation are equal, 
that is, eq. (jl]) holds for the coherent states p = \z){z\. To verify this fact, we need some further 
notations. 

For (g,p) and {u,v) in we denote z = {q + ip)/^/2, w = {u + iv)/\^, and 

a := gcos^ + psin6' = V2Re(ze~*^). 
If w = we~^^ = {u + iv) I \/2 it follows that u = u cos 6* + t> sin 6^ and 
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-a = V2Re{{w-z)e-''^). 



Then 
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It thus suffices to show that 
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Now 

= 0Fs/2fc22fe (see [H 7.374(6)]) 

V ' 

= i/2fe since y=V2Re ((u;-2)e-'«) =2-^/2 ((«;-2)e-'» + (to-2)e*'') 

^ ^(-1)^ (2fc)! - (-l^-.p)^ 

^ (2k)\ ■ A;! A;! ' ' ^ A;! ^ ' 

k=0 ^ ' k=0 

showing that eq. ([H]) is, indeed, vahd. 

The construction of the Husimi POM G|o) from the homodyne detection observable Eht by the 
generahzed Markov kernel ([7]), as given in ([S]), is based on a direct computation. Due to the 
informational completeness of the phase space observable G|o) one would expect that a similar 
reverse construction can also be given. In the appendix we demonstrate, however, that no 
simple inverse construction can be given. 

4. Discussion 

The measurements of the homodyne detection observable Eht and the phase space observable 
G|o) constitute the basic measurements in the theory of (continous variable) quantum tomog- 
raphy. Due to their informational completeness their measurement outcome statistics p^^^ and 
p^i") both separate states p. Together with the appropriate inversion formulas the state p can, 
indeed, be deduced both from p^^t from . Therefore, the experimental implemen- 

tations of measurements of Eht and G|o) as well as the inversion formulas p^^^ ^ p and p^i"^) t-^ p 
are most important, and they have widely been analysed in the literature. In spite of that we 
find it appropriate to add a few comments on these well-known methods and on the relevance 
of the kernel constructed above. 



The balanced homodyne detection is well-known and much used techique to measure a 
quandarture observable Qg, where 9 is the phase of the local oscillator, in coherent state \z), 
z = \z\e^^, and \z\ » 1. A rigorous quantum mechanical proof that the balanced homodyne 
detection observable tends, in the high amplitude limit \z\ oo, to the quadrature observable 
Qg has recently been worked out in but see also [18j and [20] • The informational complete- 
ness of the quadrature observables Qg,9 G [0,7r], is equally well-known; recent explicite proofs 
are given, for instance, in [3lll0]. The statistical sampling of various Qg is, again, well-known, 
see e.g. p!5lfT6]. though the combination of them under a single observable Eht is a recent result 
[1]. To our knowledge the only known measurement implementation of Eht, however, is the high 
amplitude balanced homodyne detection together with a random change of the phase 6. 

Phase space observables Gk, and, in particular, the Husimi POM G|o), are other important 
examples of informationally complete observables. Indeed, Gk is informationally complete if 
and only if the generating operator K is such that tr [VF^pfC] 7^ for (almost all) {q,p) G M^. 
This results is due to [2], see [I2] for some further technical elaborations. The Arthurs-Kelly 
model as well as the 8-port homodyne detection scheme are, again, well-known measurement 
models for phase space observables Gk, the latter one being also experimentally feasible. In 
fact, there is a recent rigorous proof that any phase space observable Gk can be obtained (as a 
high-amplitude limit of the local oscillator) of the 8-port homodyne detection observable, with 
the generating operator K defined by the state of the (single-mode) parametric field [9]. 

The generalized Markov kernel ([6]) allows one to obtain the Husimi distribution p^i") of any 
state p from the homodyne detection distribution p^'^' by the formula (jlj), that is, the phase 
space observable G|o) from the homodyne detection observable Eht by ([H]). This means, in 
particular, that any measurement of G|o) can, in principle, be reduced to a measurement of Eht. 



Appendix: the reverse problem 



To solve the inverse problem, we try to find a function NJ^{u,v) such that 





\w)(w\dqdp/(2Tr) 



where w 



(u + iv) I v^, and where the order of integration can be changed (at least for compact 



sets X C M). 

As before, we denote z = (q+ip) / V^, and z = ze~^^ = (q+ip) / so that q = q cos9+p sin 9, 
p = —qsm9 + p cos 9, q = q cos 9 — psm9, p = qsm9 + p cos 9, and dqdp = dqdp. 



Operate on both sides of by {—z\ ■ ■ - Iz). Now 

(-^|Eht(e X X)\z) = ^ [ {-e-''z\Q{X)\e~''z)de 

and 

{-e-''^z\Qo{X)\e-''z) = {-S\Qo{X)\S) = ^ [ e'''^-i-^? /^-ix+nf /2^^ ^ J_ f ^2vx-x^-e^^_ 
On the other hand, 

{-z\G\o){Z)\z) = ^ [[ {-z\w){w\z)dudv = ^ [[ e-l"l'-l"'l'+2*^'°("^)dudt;. 
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Thus, one must solve 



1 ^2px-.^-e ^ II Ar|'^.-(^,^)e-l-l'-l-l'+2^M^-)rfMc/w 
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or, equivalently, (note that 2ilm{zw) = i{pu — qv)) 



e 

TT 



If one makes assumptions that the right hand side of the above equation (a function of q and 

i-e.x 



p) and the function {u,v) N^^^{u,v)e integrable, one can take take the double 



Fourier transform: one gets 



If the last integral would exist, then one could change integration variables to get 

J2px-x^-q'^ + \z\^ _ p-W'^o^l'" rinrlr, — II „i2px-x'^-q^ + \z\^-i{qsme+pcose)u+i{qcose-psme)v 



Q^2px-x"-q- + \z\- _ ^-tpu^tqv^^^p _ II ^tZpx-x'-q-+\z\--i(qsmtl+pcosO)u+i(qcosO-psmtl)v ^~^p 

where the argument of the exponential function of the integrand is 

i2px — — (f + |5p — 2(g sin 6 + pcos9)u + i{q cos 9 — psin 9)v 

= -'i^ + + ^Q'(~'wsin6' + f cos6') + i'p(2x — ucos9 — t>sin6') — x^. 

Certainly the above integral does not exist even as a tempered distribution due to the |p^-term; 
for example, if u = v = x = 0, then one would have 

which clearly does not exist. To conlude, there is no function nI'^^{u,v) with the above prop- 
erties which would allow one to write Eht(0 x X) = JJ^^ Jq N^'^^ {u, v)d9dx dG\o)\{w). 
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